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SUMMARY 

A  method  is  developed  consistent  with  the  assumptions  of  small 
perturbation  theory,  which  provides  a  means  of  determining  the  down- 
wash  behind  a  wing  in  supersonic  flow  for  a  known  load  distribution. 
The  analysis  is  based  upon  the  use  of  supersonic  doublets  which  are 
distributed  over  the  plan  form  and  wake  of  the  wing  in  a  manner 
determined  from  the  wing  loading. 

The  first  application  of  the  method  proves  the  equivalence  in 
subsonic  and  supersonic  flow  of  the  downwash  at  infinity  correspond¬ 
ing  to  a  given  load  distribution.  The  principal  application  in  this 
report  is  concerned  with  the  downwash  behind  a  triangular  wing  with 
leading  edges  swept  back  of  the  Mach  cone  from  the  vertex.  A  complete 
given  along  the  center  line  of  the  wake  and  an  approxima¬ 
tion  provided  for  points  in  the  vicinity  of  this  line. 


INTRODUCTION 

The  linearization  of  the  partial  differential  equation  satisfied 
by  the  velocity  potential  for  compressible  flow  yields,  for  subsonic 
flight  speeds,  an  elliptic-type  equation  which  is  reducible  by  means 
of  an  elementary  transformation  to  the  basic  equation  in  Incompressible 
flow.  As  a  consequence  of  this  result,  wing  theory  in  the  subsonic 
realm  employs  the  same  concepts  and  same  types  of  analyses  that 
belong  to  classical  incompressible  theory.  At  supersonic  speeds,  the 
differential  equation  for  the  velocity  potential  is  hyperbolic  In  type 
and  for  wing  theory  is  equivalent  mathematically  to  the  two-dimensional 
wave  equation  of  physics.  In  spite  of  the  different  character  of  the 
basic  differential  equation  in  the  two  flight  regimes,  certain  formal 
equivalencies  can  be  set  up  which  are  intuitively  useful  in  the  solu¬ 
tion  of  specific  problems.  In  particular,  the  velocity  potentials 
of  a  three-dimensional  source  and  of  a  doublet  each  have  analogous 


forms  in  the  two  cases.  The  solution  of  different  boundary-value 
problems  encountered  in  wing  theory  has  been  discussed  in  reference 
1*  and  it  has  been  shown  how  suitable  distributions  of  sources  and 
doublets  may  be  used  to  determine  the  flow  potential  associated 
with  a  given  lifting  or  nonlifting  wing. 

The  calculation  of  downwash  behind  a  wing,  for  incompressible 
flow,  relies  almost  exclusively  on  the  use  of  Prandtl's  lifting-line 
theory  which  is,  in  turn,  developed  from  the  concept  of  a  single 
horseshoe  vortex.  The  conventional  approach  to  the  general  down- 
wash  problem  is  to  determine,  first,  the  induced  field  of  the  simple 
horseshoe  vortex  by  means  of  the  Biot-Bavart  law  and,  then,  from  a 
knowledge  of  the  spanwise  distribution  of  loading  over  the  wing,  to 
calculate  finally  the  induced  field  produced  by  a  vortex  sheet 
composed  of  superimposed  vortices  of  varying  span. 

When  downwash  calculations  are  to  be  extended  to  the  case  of 
supersonic  wings,  it  appears  at  first  that  the  use  of  vortex  sheets 
is  completely  inadmissible  since  no  equivalent  to  the  Biot-Savart 
law  exists.  It  is,  in  fact,  true  that  the  horseshoe  vortex  no 
longer  plays  the  outstanding  role  it  has  at  low  speeds.  However, 
when  a  more  detailed  investigation  is  made  of  the  underlying  analysis, 
it  becomes  apparent  that  vortex  theory  and  the  Biot-Savart  law  can 
be  developed  from  the  initial  use  of  a  constant  distribution  of 
doublets  over  a  given  surface.  (E.g.,  see  references  2  and  3.)  These 
doublets  produce  a  discontinuity  in  the  velocity  potential  at  the 
surface,  and,  for  incompressible  theory,  the  curve  which  bounds  the 
surface  can  be  identified  with  a  vortex  curve  possessing  circulation. 
The  proof  of  the  Biot-Savart  law  and  the  introduction  of  vortex  sheets 
are  direct  consequences  of  these  basic  ideas. 

Since,  as  was  shown  in  reference  1,  supersonic  boundary— value 
problems  involving  sources,  sinks,  and  doublets  can  be  solved  in  a 
manner  analogous  to  that  used  in  low— speed  theory,  a  method  is 
therefore  provided  for  an  attack  on  the  downwash  problem  for 
supersonic  plan  forms  through  the  use  of  doublet  distributions. 

The  present  report  has  two  principal  aims:  First,  to  outline 
in  some  detail  the  theoretical  approach  to  the  determination  of  the 
velocity  potential  of  the  flow  field  associated  with  a  supersonic 
lifting  surface  and  the  subsequent  calculation  of  the  downwashj  and, 
second,  to  apply  the  theory  to  the  case  of  a  triangular  wing  swept 
back  of  the  Mach  cone  and  to  present  the  results  of  the  complete 
calculations  along  the  center  line  and  in  the  plane  of  the  wing. 
Downwash  immediately  back  of  the  trailing  edge  and  at  an  infinite 
distance  behind  the  wing  will  also  be  derived  and  is  in  agreement 
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V  with  the  previously  published  results  of  P.A.  Lagerstrom  (reference 
4).  The  final  portion  of  the  analysis  will  treat  the  variation  of 
downwash  near  the  center  line  of  the  triangular  wing. 

In  the  theoretical  portion  of  the  report,  the  boundary-value 
problem  will  be  introduced  and  the  solutions,  obtained  from  Grreenfs 
theorem,  will  be  given  for  low— speed  and  supersonic  flow.  In  the 
section  of  the  report  devoted  to  applications,  the  theory  will  be 
used  first  to  evaluate  the  potential  function  at  an  infinite 
distance  downstream  from  a  lifting  wing.  Since  the  mathematical 
problems  arising  in  the  physically  obvious  case  of  the  unswept  wing 
of  infinite  span  correspond  closely  to  those  for  the  more  general 
case,  the  theory  is  next  applied  to  this  case.  From  this  applica¬ 
tion,  a  general  procedure  is  developed  for  treating  wings  with 
supersonic  trailing  edges.  The  final  application  of  the  report 
will  be  devoted  to  the  triangular  wing.  In  all  of  these  applications, 
it  will  be  seen  that  the  analytic  expressions  which  have  been  obtained 
in  supersonic  theory  for  the  load  distributions  over  certain  plan 
forms  afford  a  means  whereby  the  chordwise  distribution  of  pressure 
may  be  introduced  into  the  analysis,  and,  therefore,  such  expedients 
as  lifting— line  theory  are  no  longer  so  essential. 

No  attempt  will  be  made  here  to  discuss  the  effect  of  airfoil 

V  thickness  on  the  downwash  distribution,  although  this  effect  is 
actually  simpler  to  treat  mathematically.  It  suffices  to  state 
that  the  entire  theory  is  postulated  on  the  assumptions  of  thin— 

*  airfoil  or  small— perturbation  theory  and  that,  consequently,  thick¬ 

ness  effects  and  lifting— plate  solutions  are  additive.  For  the 
results  that  are  given  in  the  plane  of  the  airfoil,  the  thickness 
effect,  which  is  necessarily  symmetrical  with  respect  to  this  plane, 
is  zero. 


LIST  OF  IMPOBTANT  SIMBOIS 
velocity  of  sound  in  the  free  stream 
wing  span 
root  chord  of  wing 

oomplete  elliptic  integrals  of  the  second  kind  with 
moduli  k,  k0,  ki,  k2,  respectively 
2aVo 
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complete  elliptic  integrals  of  the  first  kind  with 
moduli  k,  k2,  respectively 

free—stream  Mach  number 


static  pressure 
Pz-Pu 

dynamic  pressure 


sp°T°2) 


1 / ( x-Xi ) 2+ ( y-y ! ) 2+ ( z-z  ! ) 2 

]/  (  x-Xjl  )  2+  ( y-y  x )  2+z  2 

1/ ( x-xx ) 2-P  2Cy-yi ) 2-p2( Z-Zi ) 2 


y (x— Xx) 2— P2(y~ yi) 2_P2z2 


u,v.,w  perturbation  velocity  components  in  the  direction 

of  the  x — j  y— ,  z-axes,  respectively 


Aua  uu"“uZ 

Y0  free—stream  velocity 


Wp 


z-component  of  velocity  induced  by  doublet  distribution 
over  plan  form 


wW 


z-component  of  velocity  induced  by  doublet  distribution 
over  wake 


x*y^z  Cartesian  coordinates  of  an  arbitrary  point 

xi,yi,zi  Cartesian  coordinates  of  source  or  doublet  position 
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angle  of  attack 
a/m2— 1 


semivertex  angle  of  triangular  wing 

Pyi 

CO0O 

P  tan  6 


Mach  angle {  arc  sin  — 

\  My 

density  in  free  stream 


perturbation  velocity  potential 

®u r®i 


sign  denoting  finite  part  of  integral  (equations 
(10)  and  (ll)  ) 


Subscripts 


conditions  on  upper  portion  of  surface 
conditions  on  lower  portion  of  surface 
conditions  at  leading  edge 
conditions  at  trailing  edge 


refers  to  wake 


refers  to  plan  form 

conditions  on  discontinuity  surface  (at  Zi=0) 
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Regions 

I, II, III  integration  regions  on  plan  form  (fig.  4) 

A,B  regions  in  wake  of  triangular  wing  (fig.  5) 


THEORY 

Boundary  Conditions 

The  proposed  problem  is  one  of  finding  the  downwash  behind  a 
flat  plate  which  supports  a  loading  consistent  with  its  angle  of 
attack  and  plan  form,,  It  will  be  assumed  throughout  the  analysis 
that  this  load  distribution  is  known.  Such  values  were  given  for 
several  plan  forms  in  reference  5  and.  further  results  can  be  found 
in  the  literature  on  supersonic  wings. 

The  load  distribution  over  the  wing  may  be  obtained  from  a  knowl¬ 
edge  of  the  differences  in  pressures  acting  on  the  lower  and  upper 
surfaces.  Moreover,  in  thin-airfoil  theory,  where  boundary  conditions 
are  given  in  the  z  =  0  plane  (i.e.,  the  plane  of  the  wing),  a  simple 
relation  exists  between  local-pressure  coefficient  and  the  streamwise 
component  of  the  perturbation  velocity.  Thus,  assuming  that  the  free—  * 

stream  direction  coincides  with  the  positive  x-axis  (fig.  l),  and 
denoting  by  u  the  x-component  of  the  perturbation  velocity,  it 

follows  that  r 


Ap 

q 


=  _  (Uu_uz ) 
q  v0 


2Au 
_ s 


(1) 


where  the  variables  are  defined  in  the  table  of  the  symbols. 
Furthermore,  from  the  definition  of  the  perturbation  velocity  poten¬ 
tial  £■ 


<!>  =  f"  udx  (2) 

JSi 

N. 

where  a  is  a  point  in  a  region  at  which  the  potential  is  zero. 
Combining  equations  (l)  and  (2),  the  jump  in  potential  in  the  plane 
z=0  can  be  determined  by  integrating  the  jump  in  the  u  induced 
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*  velocity  or,  what  amounts  to  the  same  things  the  jump  in  load  coeffi— 

c  ient .  Thus , 


r  vn  rx  / Ap  \ 

A4g  =  /  AugCbti  =  —  /  (  —  )  dJCi 

^  L.E.  2  J  L.E.  \  q  / 


where  the  integration  extends  from  the  leading  edge  to  the  point  x 
and  A0S  represents  the  jump  in  0  in  the  xy-plane.  Since  load 
coefficient  Ap/q  must  he  zero  off  the  wing  and  since  u  is  an  odd 
function  in  z,  the  value  of  u  must  he  zero  for  all  points  off 
the  wing  in  the  xy— plane.  It  follows  that  A0S  remains  constant 
at  a  given  span  station  for  all  values  of  x  heyond  the  trailing— 
edge  position. 

Figure  1  indicates  an  arbitrary  lifting  surface  in  the  z=0 
plane  together  with  the  distribution  of  A0S  for  given  constant 
values  of  y  and  x.  In  both  subsonic  and  supersonic  theory,  the 
wing  together  with  the  semi— inf inite  strip  extending  downstream  of 
the  wing  form  a  discontinuity  surface  for  the  velocity  potential, 
while  A$s  is  equal  to  0  throughout  the  remaining  portion  of  the 
*  xy— plane.  These  conditions,  together  with  the  fact  that  the  vertical 

induced  velocity  w  is  a  continuous  function  at  z=0,  are  sufficient 
to  determine  0  throughout  space.  The  values  of  u,  v,  and  w  can  then 
%  be  found  from  the  corresponding  partial  derivatives  of  0  with  respect 

to  x,  y,  arid  z.  The  attention  in  the  present  report  is  centered  on 
w,  the  downwash  function. 


Solution  to  Boundary— Value  Problem 

In  reference  1,  the  solutions  for  b oundary— value  problems  of 
the  type  under  consideration  were  given  for  both  incompressible  and 
supersonic  theory.  The  basic  differential  equations  satisfied  by 
the  perturbation  velocity  potential  are,  for  the  two  cases,  respectively. 


=  0 
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&c2 
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and. 


2  82iJ>  d2$  d2<I> 

8x2  8y2  8z2 


15) 


Incompressible  theory.—  For  boundary  conditions  prescribed  in 
the  z  =  0  plane;  the  solution  of  equation  (h)  is 


• (w)  ■-&//[£(!£-£*) 


taidyi  (6) 


where 


=  V^c-x  1 )  2+(y-y  1 )  2+(z-Zi ) 2 


and  t  is  the  area  for  which  the  integrand  does  not  vanish.  The  * 

terms  are  e<5ua,l  to  the  velocity  potential 

at  x,y,z  of  a  unit  source  and  doublet  situated  at  the  point 
xijyijO*  The  remaining  terms  in  the  integrand,  which  determine  the 
distribution  of  source  and  doublet  strengths,  must  be  found  from 
known  boundary  conditions.  If  a  lifting  surface  fixes  the  boundary 
conditions,  induced  vertical  velocities  on  the  upper  and  lower 
faces  of  the  surface  are  equal  so  that 

d*u  _ 

dzi  dzi 


and 


.  <  W>  -  £-/t  '-)a  ^  17) 

Equations  (6)  and  (7)  are  well  known  in  potential  theory  (refer¬ 
ence  3>  P.  60),  but  the  derivation  usually  employs  the  assumption 
that  the  value  of  Q  is  zero  at  all  points  infinitely  distant  from 
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the  wing.  This  assumption  cannot,  of  course,  he  made  in  aerodynamic 
applications  where  the  discontinuity  surface  r  extends  to  x  =  °o, 
as  in  the  case  of  a  lifting  wing  or  lifting  line  with  trailing 
vortices.  These  latter  problems,  with  which  this  report  is  directly 
concerned,  are  of  such  a  nature,  however,  that  the  induced  effects 
at  an  infinite  distance  are  confined  to  the  plane  x  =  <».  An  investi¬ 
gation  of  the  derivation  of  equation  (6)  reveals  that  the  conditions 
imposed  on  <1,  in  general,  can  he  relaxed  sufficiently  to  permit  a 
discontinuity  in  a  strip  of  finite  width  along  the  entire  extent  of 
the  x— axis.  The  mathematical  details  of  the  derivation  will  not  he 
given  here  hut  a  statement  of  the  restrictions  on  $  at  infinity  is 

A 

worthwhile.  Thus,  denoting  by  —  the  directional  derivative  of  Q 

bn 

taken  normal  to  a  prescribed  surface,  the  following  conditions 
apply: 


b® 

1.  The  functions  0  and  —  are  zero  at  all  points 

on 

having  radius  vectors  which  make  finite  (nonzero)  angles 
with  the  positive  x-axis,  the  points  lying  on  a  spherical 
surface  of  infinite  radius  with  center  at  the  wing.  (This 
preserves  the  usual  potential  theory  assumptions  except 
over  the  portion  of  the  spherical  surface  which  forms  the 
plane  x  =  «.) 


2.  The  values  0  and  —  are  bounded  at  all  points 

ox 

infinitely  distant  from  the  lifting  surface  and  at  a  non¬ 
infinite  distance  from  the  positive  x-axis.  (This  condition 

b® 

places  restrictions  on  the  values  of  $  and  —  in  the 
plane  x  =  «>.)  dx 


Conditions  (l)  and  (2)  are  satisfied  for  a  lifting  surface  of 
finite  span,  and  equation  (7)  is  consequently  applicable  directly 
to  the  determination  of  the  velocity  potential.  As  an  application 
of  the  equation,  suppose  a  sheet  of  horseshoe  vortices  is  situated 
as  in  figure  2  with  bound  vortices  placed  on  the  y-axis,  trailing 
vortices  extending  parallel  to  the  positive  x-axis,  and  has  a  span- 
wise  distribution  of  circulation  A$  symmetrical  to  the  xz— plane 

and  defined  for  —  —  <  y  <  Then  the  velocity  potential  corre¬ 
sponding  to  this  vortex  sheet  is  given  by  the  expression 


$U,y,z) 


A<Dsdyi 


dxi 


[ ( x-x! ) 2+ (y-yi ) 2+z  ]3/2 


o 


(8) 
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When  =  constant,  a  single  horseshoe  vortex  results.  * 

Supersonic  theory.—  For  supersonic  boundary-value  problems 
associated  with  plan  forma  as  Indicated  in  figure  1(a),  where  the 
known  conditions  are  given  in  the  z  =  0  plans,  the  general  solution 
of  equation  (5)  is  given  in  reference  (l)  in  the  form 


®  (x,y,z)  =  -  i- 


where 

rc  = 


and  the  subscript  s  on  the  parentheses  indicates  that  the  function 
is  to  be  evaluated  at  z\  =  0.  The  region  t  is  that  portion  of  the 
xiyi~ plane  bounded  by  the  leading  edge  of  the  wing,  the  lines  parallel 
to  the  x-axis  and  stemming  from  the  lateral  tips  of  the  wing,  and  the 
trace  in  the  Zi  =  0  plane  of  the  Mach  forecone  with  vertex  at  the 

point  x,y, z.  The  sign  |  is  to  be  read  "finite  part  of" 

and  was  introduced  by  Hadamard  (reference  6)  as  a  manipulative 
technique  with  the  property  that 


For  purposes  of  calculations,  this  was  modified  in  reference  1  to 


the  asterisk  indicating  that  no  upper  limit  is  to  be  substituted 
into  the  indefinite  integral,  the  latter  being  determined  as 

F  (x)  +  C 


V 
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4  where 


C  =  lim 
x  x0 


-  F(x)  ~ 

j 


1_ 

2n 


Equation  (9)  is  the  direct  analogue  of  equation  (6) 
l\  ,  1  (  5  1 

2n  \bzx  r, 


The  terms 


and  l  'g “  ~  )  are  elual  to  the  "velocity  potential  at 

x,y,z  of  a  unit  supersonic  source  and  doublet  situated  at  the  point 
while  the  remaining  terms  in  the  integrand  determine  the 
distribution  of  source  and  doublet  strength  and  are  determined  by  the 
known  boundary  conditions. 


When  the  potential  function  associated  with  a  lifting  surface  is 
to  be  evaluated. 


^U_*h 

Szx  dzx 


and  equation  (9)  reduces  to  the  form 


A$S  (dzT  4X1471  (12) 

s 

In  application,  the  region  of  integration  in  equations  (7)  and  (12) 
can  be  divided  into’  areas  occupied,  respectively,  by  the  plan  form 
and  the  wake  region.  Thus,  for  equation  (12), 


<J  (x,y,z) 


-zpg  P  P  _ A3>s  dxjdy! _ 

2ji  J  J  [(x- xx)2— P2(y-yi)2— |32z2]3/2 

plan  form 


z32 

rr 

A$s  dxxdyj. 

2n 

J  ^ake 

[(x-Xx )  2-32(  y-yi )  2-P2z  2  ] 3/2 

I 


% 


(13) 
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Equation  (13)  presents  a  formal  solution  for  the  calculation 
of  velocity  potential  and,  subsequently,  downwash  for  a  given  surface 
in  terms  of  A0S.  Since  A$s  was  related  directly  to  load  distri¬ 
bution  in  equation  (3),  it  is  apparent  that  the  various  known 
solutions  to  lifting— surface  problems  are  directly  applicable.  The 
fact  that  supersonic  theory  permits  the  determination  of  load 
distribution  in  closed  analytic  form  for  many  simple  plan  forms 
provides  a  distinct  advantage  that  is  lacking  in  subsonic  theory 
wherein  virtually  all  known  results  are  available  only  in  numerical 
form.  Thus,  theoretical  analysis  of  problems  involving  supersonic 
flight  speeds  can  be  carried  further  before  recourse  to  numerical 
methods  is  necessary. 


APPLICATIONS 

Value  of  Potential  Function  at  x  =  00 

It  is  possible  to  show,  from  equations  (7)  and  (12),  that  the 
potential  functions  corresponding  to  a  wing  with  fixed  load  distri¬ 
bution  are  idential  at  x  =  00  for  incompressible  and  supersonic 
flow.  Assuming  A$s  known,  the  values  of  ®  (x,y,z)  for  the  two 
cases  are  given,  respectively,  by  the  equations 


z  r  p  _ A<I>3  dx1dy1 

kitj  J  [(x-xi)2+(y-yi)2+2?]ys 

-plan  form 


z _  P  r  _ A3>s  dXidyx 

+  WJ  [(x-Xl)My-yi)2+z2F2 


wake 
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4  and 

®(x,y,z) 

32Z 
2rt 

Since,  however,  A$s  is  finite,  it  follows  immediately  that  for 
fixed  values  of  y  and  z  the  integrals  over  the  plan  form  in  both 
equations  approach  zero  as  x  increases  indefinitely.  Thus, 
denoting  by  the  value  of  xx  at  the  trailing  edge  of  the 

wing,  the  potential  functions  at  x  =  00  are  given  by  the  expressions 


(14a) 


and  for  the  supersonic  case 
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where  the  symmetry  of  the  problem  with  respect  to  the  yi-axis  has  been  used.  Computing  the 
finite  part  of  the  integral 
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For  a  point  'behind,  the  trai ling-edge  wave  (fig.  3("b))j  the  two  quantities 
determined  in  a  similar  manner.  Thus, 
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These  results  are  the  familiar  equations  associated  with  two-dimensional  supersonic  flat- 
plate  theory. 
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* 


The  point  of  principal  interest  in  this  development  is  the  i 

jump  in  the  induced  vertical  velocity  w  in  passing  from  a  point 
just  ahead  of  the  trailing-edge  wave  to  a  point  immediately  behind 
the  wave.  A  study  of  equations  (17)  and  (19)  shows  that  this  jump 
is  the  result  of  the  discontinuity  in  the  contribution  to  the  down- 
wash  of  the  term  containing  the  single  integral.  Ahead  of  the 
trailing-edge  wave,  this  term  yielded  the  result  that 


w  =  —  V  a 
0 

whereas  behind  the  wave,  the  contribution  of  the  term  to  w  was 
zero. 


The  method  of  attack  used  in  the  study  of  the  unswept  wing  can 
be  generalized  to  apply  to  arbitrary  plan  forms.  A  discussion  of 
this  case  follows. 

Arbitrary  plan  forms.—  As  will  be  shown  later,  for  any  plan 
form  with  supersonic  trailing  edge,  the  jump  in  the  value  of  w 
in  the  plane  of  the  wing  at  the  trailing  edge  can  be  calculated 
directly  by  means  of  simple  momentum  methods.  At  this  point,  how¬ 
ever,  it  is  of  more  interest  to  consider  in  a  genral  manner  the  „ 

nature  of  the  integrations  involved  when  the  point  x,y,z  is  either 
ahead  of  or  behind  the  trailing— edge  wave.  Figures  4(a)  and  4(b) 
show  a  plan  form  with  a  straight  trailing  edge  with  areas  of  integration  * 

indicated  for  the  point  P  in  each  of  the  two  positions.  (The 
straight  trailing  edge  is  not  a  necessary  restriction  and  is  only 
introduced  for  convenience  of  notation.)  The  regions  of  integration 
are  divided  under  the  assumption  that  the  first  integration  on  the 
plan  form  in  equation  (13)  will  be  made  with  respect  to  y1#  When 
the  point  P  is  ahead  of  the  trailing-edge  wave,  therefore,  the 
contribution  of  the  wake  is  zero  and  the  integration  over  the  plan 
form  is  made  to  conform  with  regions  I  and  II.  When  the  point  is 
behind  the  trailing-edge  wave,  three  integrals  are  evaluated 
corresponding  to  regions  I,  II,  and  III.  In  the  case  of  the  infinite 
aspect  ratio,  unswept  wing  region  I  was,  of  course,  nonexistent  and, 
in  general,  no  essential  difficulty  in  regard  to  the  limits  of 
integration  is  introduced  by  this  region  regardless  of  where  P  is 
situated.  In  region  II,  however,  the  problem  must  be  treated  in 
more  detail. 

Consider  first  the  case  when  P  is  ahead  of  the  wave  and  denote 
by  ^Ila  the  contribution  of  one  side  of  region  II  to  the  total 
potential.  Then 
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Similarly,  when  P  is  behind  the  wave  and  the  same  subscript 
notation  is  used  to  refer  to  one  side  of  region  II,  the  value  of 

®IIa  ^ 


— 02z 

®  II a  =  lim  - 

o  2  it 


r  1  ^>1  (xxJyi)dy1 _ 

y  +  e  [  ( x-Xi )  2-02  ( y-yi  )2-02z2]3/2 


(21b) 


where  Yi  is  as  defined  above. 

The  contributions  of  the  other  side  of  region  II  to  the  potential 
will  not  be  considered  separately  as  the  behavior  is  identical.  When 
P  lies  ahead  of  the  wave,  e  appears  in  the  limits  for  integration 
with  respect  to  both  jx  an&  x1#  This  corresponds  to  the  situation 
in  equations  (l 6)  and  (17)  and,  as  for  that  problem,  the  limiting 
process  is  carried  out  after  the  integration  is  completed.  When  P 
is  behind  the  wave,  it  is  not  necessary  to  defer  the  limiting  process, 
since 
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But  since 


- - -  ^  M 

[  ( x— x  1 ) 2— P  2  ( y— y  1 ) 2— P  2z  2  ] 3  /2  = 


for  c0  ^  xx  ^  x0  (i.e.,  J 
interval  of  the  first  integra: 


.  is  hounded,  for  all  values  of 
ion);  and  further,  since 


Xi  in  the 


J+e 

lim  M  /  dyx  =  lim  Me  =  0 

£  O  Jy  £  ->  0 

therefore,  for  P  situated  behind  the  trailing— edge  wave,  the 
contribution  of  region  Ila  is  given  by 


_ (x1,y1)dy1 _ 

[  ( x-xi ) 2— P2(  y-y  1  )2~P2z2]  3/£ 


(22) 


¥ 


Equation  (22)  will  “be  applied  directly  in  the  determination  of 
downwash  behind  the  triangular  wing.  The  significance  of  the  result 
is  that,,  when  the  point  P  at  x,y, z  is  behind  the  Mach  wave  from 
a  supersonic  trailing  edge,  the  limiting  process  associated  with 
region  II  need  not  be  considered.  When  P  is  ahead  of  the  Mach  wave, 
the  term  e  must  be  retained  in  the  analysis  and  the  limiting  process 
used.  As  was  previously  noted,  the  general  analysis  developed  in 
this  report  places  no  restriction  on  the  orientation  of  the  trailing 
edge;  however,  it  should  be  pointed  out  that  region  II  exists  only 
for  the  case  in  which  the  trailing  edge  is  supersonic.  Therefore, 
the  jump  in  downwash,  obtained  from  the  integration  over  region  II, 
is  associated  only  with  supersonic  trailing  edges;  whereas  both  the 
downwash  and  loading  are  continuous  across  a  subsonic  trailing  edge. 


Triangular  Wing 

Consider  a  triangular  wing  (fig.  5)  with  leading  edges  swept 
back  of  the  Mach  cone  from  the  vertex.  The  loading  over  the  wing 
is  known  to  be  (references  7  and  5) 


if 
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X 


Ap  40o2(xxi 

q  E0P  V0o2Xi2-P2yi2 


(23) 


where  E0  is  the  complete  elliptic  integral  of  the  second  kind  with 

modulus  k0  =  (/l—Q02  and  Q0  =  (3  tan  6,  5  "being  the  semivertex 

angle  of  the  triangle.  From  equation  (3) 

A$a  =  H  't/e02x12-p2j12  (24) 

where 


2aVo 

E0P 


(25) 


Sufficient  information  is  now  at  hand  to  permit  the  use  of  equa¬ 
tion  (13).  Setting,  for  convenience,  $=  Op  +  Oy  the  velocity 
potential  at  x,y,z  is  given  "by  the  sum  of  the  two  expressions 


Op  =  — 


zHp 

2* 


yl'go2xig-P2yig  foe!  Ay! 

plan  form  [  (x-Xi)2-p2(y-yi)  2-p2z2]3/2 


(26) 


%  =  “ 


zHP 

2rt 


/Le 


/e02c02-P2yi2  dxi  dyx 
[  (x-xx )  2-p2(y-yx )  2-p2z2] 3/; 


(27) 


Equation  (26)  represents  the  contribution  to  the  velocity  potential 
furnished  "by  the  doublets  distributed  over  the  plan  form, while 
equation  (27)  represents  the  contribution  furnished  by  the  doublets 
within  the  wake.  The  latter  equation  is  the  mathematical  equivalent 
in  supersonic  flow  of  the  subsonic  velocity  potential  of  a  sheet  of 
horseshoe  vortices  corresponding  to  an  elliptic  span  load  distri¬ 
bution.  Equations  (l4a)  and  (l4b)  showed  that  the  expression  for 
at  x  -  00  is  identical  to  the  velocity  potential  of  the 
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subsonic  vortex  sheet.  However,  in  the  vicinity  of  the  x  =  c0  * 

line,  the  behavior  is  entirely  different. 

In  the  present  report,  equations  (26)  and  (27)  will  be  applied 
only  to  the  determination  of  downwash  along  the  center  line  in  the 
wake  of  the  airfoil.  Further  extensions  are  limited  only  by 
difficulty  of  integration.  It  will  appear  in  the  development  that 
recourse  to  numerical  methods  is  probably  necessary  in  general. 

Downwash  induced  by  doublets  in  the  wake.-  Setting  y  =  0  in 
equation  (27)  and  integrating  with  respect  to  Xi,  it  follows  that 


zH(x-Cq)  V^o2Cq2— P2yi^  dyx _ 

2 It  — Yi  (yi2+z2)  fj ( x— c  o )  2— P 2y  1  2  P 2z  2 


(28) 


where  the  limits  on  the  integral  are  not  yet  specif  ied,  since  they 
differ  in  the  regions  A  and  B  indicated  in  figure  5*  In  either 
case,  however,  the  limits  are  seen  to  be  the  roots  of  one  of  the  two 
radicals  in  the  integrand. 

In  order  to  derive  an  expression  for  downwash  in  the  plane  of 
the  airfoil,  it  is  convenient  to  express  equation  (28)  in  a  different 
form.  Integrating  by  parts 


$W  =  Ve02Cog-p2yi 
H  2it 


arc  tan 


yi(x-c0) 


z  /(x-c0)2H32yi2-P2z: 


-Ya 


yi 

Ve02c02-P2yi2 


arc  tan 


yiU-Cp) _ 

z  hj ( x  c o )  ~  {3 2y  1  2H32z2 


^  w 

This  form  of  the  equation  is  easier  to  treat  when  lim  - —  is  to 

z->  ooz 

be  considered.  In  both  regions  A  and  B,  it  can  be  shown  that  wy, 
the  contribution  to  the  downwash  made  by  the  doublets  in  the  wake, 
is  given  by  the  expression 

if 
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H32 

2jt  ( x— c  o ) 


0o2co2-P2yi 


d-ya 


(29) 


In  region  B,  Yi  = 


60Ci 

"T 


Using 


an  even  function  and  introducing  the 


the  fact  that  the  integrand  is 
notation 


Setting  kx  =  — — 
x0~l 

expression  for  Wy 


and  noting  that  ki  <  1  in  the  region,  the 
in  region  B  is 


(31) 


where  Ei  is  the  complete  elliptic  integral  of  the  second  kind  and 
the  subscript  denotes  that  the  modulus  is  ki . 

In  region  A 


Yi 


P 


By  means  of  the  transformations 

x 


•a  =  pyr. 

x-^o 


X 


2k 


MCA  TN  No.  1620 


the  integral  is  transformed  to  the  form 


# 


HP 


1 


dr| 


Setting  k2  =  -  = 

#o 

introducing  Jacobian 
T|  =  sn  u  ,  we  have 


1 

— ,  noting  that  k2  <  1  in  region  A,  and 
ki 

elliptic  functions  in  the  transformation 


Wijrt  n  K2  (l— k22)K2-E2 

— —  =  -  k2  /  cn2ud.u  =  - - -  (32) 

HP  JQ  k2 


where  K2  and  E2  are  complete  elliptic  integrals  of  the  first  and 
second  kind  with  modulus  k2. 


The  values  of  the  downwash  given  by  equations  (31)  and  (32)  can 
be  presented  in  terms  of  w0,  the  downwash  on  the  wing,  where  from 

— w0 

equation  (25)  and  from  the  boundary  condition,  a  =  - — , 

^  o 


w0tc 

W 


(33) 


Thus, 

Wy.  E2—(l—k22)K2 

—  =  2  -  (  region  A) 

w0  jtEok2 


nr  _  2Ei 
Vo  ”  5®^ 


(  region  B) 


ww 

Figure  6  shows  this  variation  in  —  plotted  as  a  function  of  xQ 

w° 

for  values  of  0O  equal  to  0.2,  0.4,  0.6,  and  0.8. 

Downwash  induced  by  doublets  on  the  plan  form.—  Along  the  center 
line,  the  contribution  to  the  potential  function  due  to  the  distribution 
of  doublets  on  the  plan  form  is  given  by  equation  (26) ,  where  the  value 
of  y  is  zero.  Just  as  in  the  case  of  the  wake  distribution,  it  is 
seen  that  the  nature  of  the  solution  changes  in  passing  from  region  A 
to  region  B  so  that  the  derivations  for  these  regions  will  be  given 
independently . 


¥ 


% 


MCA  TN  No.  1620 


25 


In  region  B,  no  singularities  occur  in  equation  (26)  so  the 
finite  part  sign  can  be  disregarded.  Since,  moreover,  the  integrand 
is  an  even  function  of  yl3  the  expression  for  velocity  potential 

D  ©  C  Oili©  3 


After  changing  variables  by  the  transformations 


equation  (3^-)  becomes 


IX 


°  *L  r1  _jni 

9°  1  Jn  (l-k2!2)1 


Substituting  £=sn  u  and  noting  that  sn  K=l,  sn  0=0  where  K  is 
the  complete  elliptic  integral  with  modulus  k,  we  have 


Op  =  — 


■co  >3 


eV  dxx  r  cd2udu  =  -  f  °  -  QHQdxx 
’’  1  o  n  Jo  */(x-Xi)2— 32Z2 


Tne  dqwnv/ash  wP  in  the  z=0  plane  can  now  be  found  by  considering 

lim  — This  pea^a  the  result 
z^  o  d  2 


where  the  modulus  reduces  to 
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0Q*1 
X— Xi 


Kewriting  the  equation  so  that  k  is  the  variable  of  integration 

and  ki  =  0  -- 

xQ-l 


(35) 


As  is  to  he  expected,  the  value  of  is  seen  to  approach 

zero  as  xQ  becomes  infinitely  large.  The  upstream  boundary  of 
region  B  lies  at  Xo=l+0o  and,  as  will  be  seen,  Wp  is  continuous 
at  this  point. 

Consider  now  region  A.  In  this  case,  the  traces  of  the  Mach 
forecone  of  the  point  (x,o,z)  cut  across  the  plan  form.  For  ease 
of  calculation,  the  area  within  the  wing  leading  and  trailing  edges 
and  the  traces  of  the  forecone  will  be  divided  into  two  parts.  Thus, 
making  use  again  of  the  symmetry  with  respect  to  the  xi-axis, 
and  also  the  results  discussed  in  connection  with  equation  (22),  it 
follows  from  equation  (26)  that 


The  reduction  of  the  first  double  integral  in  equation  (36)  can  be 
carried  out  by  methods  exactly  equivalent  to  those  used  in  reducing 
equation  (3^)  to  (35)*  As  a  result  of  this  calculation,  the  contribu¬ 
tion  made  to  Wp  in  the  plane  z=0  by  the  first  double  integral  is 
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e 


_  h pr1  K— 1 

*Jn  k+eQ 


dk 


(37) 


In  the  second  double  Integral  of  equation  (36)  the  following 
substitution  is  made 


I  = 


Pyi 


y^x-Xi)2-^ 


k  = 


J ( x— Xi ) 


2-£2z2 


0OX1 


Then^  if 


I  = 


J/T 


-!Ll)^Pa22 


2xi2-p2yi2 


[ ( x— xx ) 2— p2y x2  P2z2] 


d-yi 


I  can  be  written 

t  -  1  r  1  ^--k2|g 

P0oXlk2J  (1-S2)3/2 

o 


The  integration  can  "best  he  completed  here  through  the  use  of  equa¬ 
tion  (l4)  where 


A(|) 


(l+|)3/2 


and  the  auxiliary  transformation  5  =  sn  u  is  used.  As  a  result  of 
this  calculation 


I  =  — ~ — x  (k-e) 

pe0xik2 


* 


so  that  the  second  double  integral  of  equation  (36)  is  equal  to 


zHp  n  °o  K-E 
Jt  eoxik2 


* 
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Using  methods  similar  to  those 
tion  (35) 1  the  contribution  to 

x^"“X 

where  k2  =  - .  Equations  (37)  and  (38)  may  now  be  combined  to 

®o 

give  the  total  downwash  in  region  A  induced  by  the  doublets  on  the 
plan  form.  Thus , 


employed  in  the  derivation  of  equa— 
w-c  is  found  to  be 


K-E 


^(1+0  Qk) 


dk 


(38) 


# 


wpir 

HP 


K-E 

k^l+Gok) 


K-E 

k+0o 


dk 


(39) 


The  integration  of  equation  (39)  requires  numerical  methods.  A 
slight  simplification  can,  however,  be  introduced  by  using  the  known 
value  of  the  downwash  w3  just  behind  the  trailing  edge.  As  will 
be  shown  in  the  next  section. 


w3* 

Hp 


=  -f  (EO-0O) 


Equation  (39)  must,  of  course,  yield  this  result  for  x0  on  the 
trailing  edge,  that  is,  for  Xq=1  and  k2=0.  Using  this  relation 
together  with  equation  (39)  it  follows  that 


Wpjr 


(40) 


Values  given  try  equations  (35)  and  (40)  are  consistent  at  the 
point  ki=k2=l,  that  is,  at  the  point  where  the  Mach  cones  from  the 
two  trailing— edge  tips  intersect  on  the  x— axis.  By  means  of  these 
equations,  the  contribution  to  the  downwash  of  the  doublets  distri¬ 
buted  over  the  plan  form  is  determined. 
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In  summary,  these  relations  are 


=  Eo-^o  _  _2_  rk2  K-E 
wQ  E0  JtE0  J  k2(l+0ok) 

o 

25  _  jl.  rkl 

wo  «E0  J0 


(region  A) 


(region  B) 


Figure  7  shows  the  results  of  the  integrations,  the  function  — 

wc 

being  plotted  as  a  function  of  Xq  for  various  values  of  dQ. 

Conditions  at  the  trailing  edge.—  The  value  of  the  vertical 
induced  velocity  immediately  ahead  of  and  behind  the  trailing-edge 
wave  must,  of  course,  be  determinable  directly  from  equations  (21 ) 
and  (22),  respectively,  by  setting  x=c0  ±  z|3.  If,  however,  the 
discussion  is  restricted  to  the  z=0  plane,  a  much  simpler  method 
exists  for  finding  w  at  these  points.  The  approach  taken  here 
follows  essentially  that  given  by  Lagerstrom  in  reference  k. 

Let  conditions  Just  ahead  of  the  trailing-edge  wave  be  denoted 
by  the  subscript  2  and  conditions  Just  behind  the  wave  by  -the 
subscript  3.  Figure  8  shows  a  section  of  a  given  wing  in  the  plane 
y=constant.  The  Mach  waves  at  the  leading  and  trailing  edge  make  the 

angle  p=arc  sin  i  with  the  z=0  plane,  and  the  wing  is  presumed 

to  be  at  angle  of  attack  a.  Assuming  the  trailing  edge  to  be 
normal  to  the  free-etream  direction,  the  variation  in  the  x-component 
in  velocity  when  passing  through  the  trailing-edge  wave  can  be  treated 
as  a  two-dimensional  problem  with  the  condition  imposed  that  U3=0 
in  the  z=0  plane .  It  is  known  that  continuity  of  flow  together  with 
balance  of  tangential  momentum  across  the  wave  lead  to  the  result 
that  the  component  of  velocity  tangential  to  the  wave  is  continuous. 
The  tangential  components  of  velocity  immediately  ahead  of  and  behind 
the  wave  are  given,  respectively,  by  the  expressions 


(vt)  =  (V0+u2)  cos  p  +  w2  sin  p 
(vt)a  =  V0  cos  p  +  w3  sin  p 
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Equating  these  relations, it  follows  that 


w3  =  w2  +  u2  cot  |i 

VoP/Ap\ 

■-roa*—(rl 


(hi) 


For  the  two  plan  forms  which  have  "been  considered,  equation  (4l) 
gives  the  following  results 

(a)  Unswept  wing  of  infinite  span 


In  this  case 


and 


4a 

0 


W3 


-  V0a 


4a 

T 


0 


(42) 


Equation  (20)  showed  that  this  result  actually  applies  to  all  points 
behind  the  Mach  wave . 

(b)  Triangular  wing 

From  equation  (23) 


/4p\  =  4e02ac0 

\  q4  EoP  l/Sp^co^-P^T^ 

from  which  it  follows  that 


W3  =  wo 


(43) 
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Setting  y i=0  and  introducing  H  from  equation  (25) 


V3* 

HP 


m 


which  is  the  expression  used  in  equation  (4o),  These  results  are 
equivalent  to  those  given  "by  Lagerstrom  in  reference  k. 

Values  of  downwash  near  center  line  of  wake*—  The  values  of 
downwash,  which  were  obtained  on  the  center  line  of  the  wake,  were 
worked  out  in  exact  detail  and  subjected  to  no  restrictions  other 
than  those  originally  imposed  by  the  use  of  the  linearized  equations 
of  flow*  It  is  possible,  moreover,  to  get  an  indication  of  the 
variation  in  the  downwash  function  in  a  portion  of  the  wake  for 
points  near  the  center  line  through  the  use  of  a  generalized  Taylor's 
expansion  in  the  vicinity  of  the  line  y=z=0*  The  next  higher  terms 
in  the  expansion  can  be  found  without  too  much  difficulty  for  the 
region  bounded  as  follows 

(a)  -|b<y<|b 


(b)  Both  y  and  z  lying  within  the  Mach  cones  from  the 
trailing-edge  tips 

The  problem  resolves  itself  into  one  of  finding  the  first 
nonvanishing  coefficients  of  y0  and  z0  in  the  series 

=  A0+A1z0+B1y0+A2z02+C2z0y0+B2y02+...  (1*5) 


The  value  of  A0  has,  of  course,  already  been  computed  and  is 
known  from  equations  (31)  and  (35)  to  be 


A0 


(46) 


0q  w 

x0-l* 


where  ki  = 
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The  coefficient  Bi  in  the  expansion  is  known  to  be  equal  to 

The  evaluation  of  this  term  follows  a  devlopment  similar 

to  that  used  for  equations  (29)  and  (34)  except  that  y  is  retained 
in  the  analysis.  The  derivation  yields  the  relation 


(x0-l)2-(y0-Tl)2 

e0z-n2  u‘ 


wit 

Hp  “  2(x0-l) 


Xo-1 


1  r » r° ....  _ 4, 

2  i  it o0  I  Uo-sJMyo-if  )3'2 


Carrying  out  the  differentiation,  with  proper  regard  for  the  singu¬ 
larity  in  the  first  integral,  it  follows  -chat  Bi  =  0.  Similarly, 


it  can  be  shown  that  C2  =  0,  while  the  coefficient  B2 
is  given  by  the  expression 


B2 


1 

itlo0oe 


ki 


itEo0o2xo2  J0 


(k+en) 


E 


1+k2 


-  E 


(1-k2)2  l^2 


dk 


(47) 


where  again  ki  =  — 2—, 

Xq—1 

In  order  to  calculate  the  variation  with  z,  it  is  necessary 
to  evaluate  Ai  = 

‘O' 


r 
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where 


—  =  p  f 

HP  J  n 


0QCO 


Jl 


y/e02c02-P2yi2 


arc  tan 


yi(x-co) 


z  /(  x-c  0 )  2-P  2yi  2-P  2z  2 


dyi 


e0xi 


/l/e02xi2-p2yi2 


[ ( x-xx ) 2-p2yi2-p2z2 ] 3/2 


dyi 


The  double  integral  contributes  nothing  to  the  coefficient,  and  the 
remaining  portion  of  the  expression  can  be  evaluated  without 
integrating  by  differentiating  twice  and  using  Cauchy1  s  integral 
theorem 


Thus  , 


r  f(yi)dyi  /  df\ 

J  7^)2=  2ni  (m\1=1z 


m 


The  coefficient  A2  will  not  be  evaluated,  since  the  first 
higher  order  term  in  z  has  been  found.  Thus,  to  the  first  order 
in  v^  and  z~ ,  the  downwash  function  —  is 

V  O  U  J  TT 

wo 


i  =  —  fix  +  fkl—  6k)  -  -52-  JL 

v0  nE0  V  j0  k+eo  /  Eoeo  lzl 


(^9) 


DISCUSSION 

The  variation  of  — ,  the  downwash  due  to  the  doublets  in  the 
wo 

wake  in  terms  of  the  downwash  on  the  triangular  wing,  was  given  in 
figure  6  as  a  function  of  x0  for  various  values  of  0O.  This 
distribution  of  downwash  along  the  center  line  of  the  wake  can  be 
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presented  as  a  single  curve  if  —  E0  is  expressed  in  terms  of 

wo 

A  graph  of  this  function  is  shown  in  figure  9. 


Figure  10  presents  the  resultant  downwash  on  the  center  line 
induced  by  the  doublets  distributed  over  the  wake  and  over  the  plan 
form.  This  is  the  complete  value  of  downwash  on  the  x0-axis  and 
shows  that,  following  the  discontinuity  in  downwash  at  the  trailing 

edge,  the  magnitude  of  builds  up  to  its  asymptotic  value  and 

w0 

achieves  an  almost  constant  value  within  approximately  one  chord 
length  aft  of  the  trailing  edge  for  values  of  0o  less  than  0.8. 

The  table  given  on  the  figure  relates  90  to  free-stream  Mach  number 
M  and  to  p  for  the  particular  case  of  a  triangular  wing  with  lead¬ 
ing  edges  swept  back  45°. 

It  has  been  shown  that  at  x  =  00  the  downwash  has  the  same  form 
for  both  incompressible  and  supersonic  flow,  provided  the  span  load 
distributions  are  equal.  An  exact  solution  for  the  downwash  on  the 
center  line  of  the  triangular  wing  has  been  developed  and  extended 
to  include  a  region  near  the  center  line  of  the  wake  by  means  of  an 
approximation  (equation  (45)).  A  measure  of  the  exactness  of  the 
approximation  may  be  obtained  by  comparing  these  results  with  the 
exact  results  for  the  incompressible  case  for  an  elliptic  span  load 
distribution,  as  obtained  from  reference  8,  page  151.  Figure  11 
shows  the  exact  values  of  for  incompressible  theory  for 

y0  zD  w° 

0  <  ~  <0.3  and  0  <  - —  <  0.3.  Also  included  are  the  linear  approxi— 

* 0  .  9°  w 

mations  m  supersonic  theory  to  —  E0  for  the  same  range  of  the 

wo 

geometric  variables.  The  span  load  distributions  for  both  wings  are 
elliptical,  but,  in  order  that  the  same  value  of  lift  is  maintained 
for  the  low—  and  high-speed  cases,  the  angle  of  attack  must  be 
modified.  This  accounts  for  the  use  of  the  factor  Eq  in  one  s(et 
of  ordinates.  It  is  apparent  from  the  comparison  that,  at  a  large 
distance  behind  the  wing,  the  approximation  of  equation  (49)  yields 
results  within  4  percent  of  the  exact  value  for  the  region  considered. 


In  application,  it  is  desirable  to  know  not  only  the  order  of 
the  variables  retained  in  the  analysis,  but  also  the  magnitudes  of 
the  neglected  terms.  In  this  respect,  equation  (49)  does  not  supply 
as  much  information  as  might  be  desired,  since  no  estimation  is 
furnished  oi  further  coefficients  in  the  series.  It  appears,  however, 
from  the  nature  of  the  agreement  obtained  in  figure  11  for  conditions 
at  a  large  distance  from  the  wing  together  with  observations  as  to 
the  asymptotic  behavior  apparent  in  figure  10,  that  equation  (49) 


r 
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furnishes  a  reasonable  estimate  of  downwash  in  region  B,  provided 
that  the  point  of  calculation  is  not  in  the  immediate  vicinity  of 
the  tip  Mach  cones.  For  values  in  region  A,  no  conclusions  can  be 
drawn  from  this  investigation  about  the  variation  of  downwash  with 
respect  to  z.  The  spanwise  variation  indicated  by  equation  (49) 
agrees  with  the  value  given  by  equation  (4-3)  at  the  trailing  edge; 
that  is,  to  the  first  order  in  y,  the  downwash  is  a  constant. 

Thus,  in  the  xy— plane,  equation  (49)  can  be  expected  to  give  a 
reasonably  close  estimate  to  the  downwash  near  the  center  line  and 
away  from  the  tip  Mach  cones  for  all  points  behind  the  trailing  edge . 
For  variation  in  the  z-direction,  however,  the  validity  of  equation  (49) 
is  restricted  to  points  in  region  B. 

The  behavior  of  the  downwash  in  the  vicinity  of  the  trailing 
edge  indicates  that  in  this  region  an  important  difference  exists 
between  lifting— line  theory  in  supersonic  and  in  subsonic  flow. 

Figure  6  shows  the  variation  of  downwash  along  the  center  line  of 
a  doublet  sheet  (lifting— line  theory  in  incompressible  flow)  with  an 
elliptic  span  load  distribution.  The  error  introduced  by  not  using 
the  chordwise  distribution  of  load  is  given  in  figure  7  and  is 
sizeable  for  about  one  chord  length  back  of  the  wing. 

The  methods  of  analysis  presented  here  were  shown  in  the  first 
application  to  lead  immediately  to  the  value  of  the  potential 
function  at  x  =  co,  and  thus  provide  a  ready  method  for  the  deter¬ 
mination  of  vortex  drag  of  a  supersonic  wing. 

A  more  detailed  study  of  downwash  behind  the  wing  will  necessarily 
involve  considerable  labor.  The  methods  given  in  the  report  are 
however,  general  and  are  directly  applicable. 
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Figure  I.  -  Sketch  showing  arbitrary  lifting  surface  together 
with  distribution  of  A<f>s,  the  jump  in  perturbation  velocity 
potential  in  the  plane  of  the  surface . 


(a) Point  P  ahead  of  trailing  edge  (b)  Point  P  behind  trailing  edg 


Figure  4-  Areas  of  integration  for  arbitrary  plan  form  with 
supersonic  trailing  edge. 


Figure  5.  -  Triangular  wing  swept  behind  Mach  cone 
showing  location  of  regions  A  and  B. 
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!0}  0O=.2 


x0  ,  Distance  in  chords 


Figure  6  -  Variation  of  the  part  of  downwash  on  x-axis 
induced  by  doublets  in  wake  with  distance  downstream 
in  chord  lengths ,  xQ  =  x/cQ.  Triangular  wing. 


Figure  7.  -  Variation  of  the  part  of  downwash  on  x-axis  induced 
by  doublets  on  the  plan  form  with  distance  downstream  in 
chord  lengths ,  xQ=  x/cQ.  Triangular  wing . 
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Figure  8  -  Sketch  of  velocity  vectors  of  the  air  before  reaching, 
on,  and  after  leaving  supersonic  airfoil. 


Figure  9  -Downwash  factor  induced  by  doublets  in  wake  ,wwE0/wot 
plotted  against  factor  representing  distance  downstream 
from  trailing  edge,  (x0- !)/0o  Triangular  wing. 
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behind  Mach  cone  with  distance  downstream  in  chord  lengths ,  x  =  x/c0 . 
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approximate 


Figure  it.  -Dorr mash  at  a  large  distance  behind  a  wing  with 
elliptical  span  loading  in  either  subsonic  or  supersonic 
flight  compared  with  first  order  expansion  determined 
from  equation  (44). 


